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1. Introduction 

The method for finding approximate solutions to the Schrodinger equation based on a matrix 
representation for the Hamiltonian is widely used in nonrelativistic quantum mechanics 1 1 1. 
Its rationale is the intuitive idea that the larger a finite basis used in particular calculations, the 
more exact are the results that may be achieved, and the "most exact" solution of the problem 
corresponds to the entire (infinite) basis set. Of course, there is now a rigorous justification 
for this physically intuitive idea within functional analysis |2|, but we would like to stress 
here that up to now only a few cases are known where one is able to manipulate a whole 
Fourier series expansion and get the exact solution to the problem. In particular, the well- 
known J-matrix method |3| (although it has recently been significantly modified |4|), uses 
a special truncation procedure for modelling a real scattering process involving interacting 
particles. Till now this method uses mainly the possibility of getting exact solutions to the 
simplest quantum mechanical problem, namely, solving the free particle Schrodinger equation 
expressed in a special basis as the eigenvalue problem for a tridiagonal (Jacobi) matrix. 

Recently |5| it has been shown that an interaction can be included in such a model 
without affecting the tridiagonal form of the Hamiltonian, a fact that can significantly enlarge 
the range of possible applications for the J-matrix method. In this paper we generalize the 
results of 1 5 1 to the case of a matrix Schrodinger equation, obtaining new exactly solvable 
multichannel potentials and, at the same time, opening the way for much wider applications 
of the J-matrix method. In particular, we will show that there exist exactly solvable potentials 
whose matrix representation has the form of a tridiagonal matrix whose non-zero entries are 
not numbers, but matrices. 

The paper is organized as follows. In the next Section we show how an eigenvalue 
problem for a Hamiltonian possessing both external and internal degrees of freedom (i.e. 
a multichannel Hamiltonian), may be reduced to the eigenvalue problem for a difference 
operator having operator-valued coefficients. We introduce a special space of operator-valued 
sequences T£ 2 , which is the counterpart of the Hilbert space I 2 . In this space we define 
an intertwining (or transformation) operator and its adjoint which relate solutions to two 
different operator-valued discrete eigenvalue problems. In Section 3 we give a solution to 
the intertwining relation, specifying the transformation operator as well as the transformed 
Hamiltonian. We also show that in the particular case of a two-dimensional space for an 
internal degree of freedom, our formulas acquire a simplest, though nontrivial, form and we 
are able to establish a one-to-one correspondence between the spaces of solutions of the initial 
and transformed equations. Section 4 is devoted to applying our method to finding a wide class 
of new exactly solvable two-channel Hamiltonians whose matrix representation has infinite 
block-tridiagonal form. In the last Section we outline a possible continuation of this work. 

2. Preliminaries 

We will show here how the usual eigenvalue problem for a Hamiltonian with a composite 
interaction may be reduced to an eigenvalue problem for a difference operator having 
operator-valued (or matrix-valued) coefficients. 

To begin with, let us consider a time-independent Hamiltonian 7io describing a 
composite system possessing both external and internal degrees of freedom. The external 
degrees of freedom are related to the motion of the system, considered as a material point, 
interacting with time-independent external fields. We shall deal only with stationary states 
and describe this motion with the help of a separable Hilbert space H which, in particular, 
may be chosen as the space L 2 (R). Internal degrees of freedom, related to the system's 
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structure (spin, colour, charm, . . . ), will be associated with an TV-dimensional linear space 
C N . A typical example is a multichannel quantum system 1 6 1 . Thus, the Hamiltonian 7i acts 
in the Hilbert space H = H <E> C N and we consider the eigenvalue problem 

H Q \V)=£\V), |*>6H. (2.1) 

Let {|n), n = 0, 1, 2, . . .} be an orthonormal basis in H and {e a , a = 1, . . ., TV} be an 
orthonormal basis in C w , so that the set {\n, a) = \n)e a = e a \n}} is a basis in the space H. 
Hence, denoting the Fourier coefficients of an element \^>) € H by Q in the basis |n, a), 
one has 

oo 

|*) = ^Vn,a|^,a) = ^*„|n) (2.2) 

where ty n — X}a=i Ai, a e a denotes a column-vector of C^. Since the right hand side of the 
equation ( 12.21 looks like a development of a vector |^) over the basis {|n)}, the quantities 
{$„} may be considered as the "coordinates" of the ket |\&). Following this approach, we can 
consider the elements | of H, as a set of elements endowed with a different structure: it is 
not a vector space anymore, since its underlying space is neither a field nor a ring (only 
the addition of this type of objects is defined, but not a multiplication by the elements of C N ); 
rather, it can be treated as a module over the linear space C , which does not correspond to 
an usual algebraic construction of a module over a ring or a field [7|. We shall denote the 
set {|^)} with this new structure by T. Although the kets {\n}} are not elements of T, we 
will call them a "basis" in the sense that V|^) G T there exists in <C N a set {^ n } such that 
d2.2i is valid, in the sense that it corresponds to a convergent Fourier series in H. Note also 
that no notion of linear dependence or independence can be introduced in T, since it is not a 
linear space, but a module. Nevertheless, one can introduce an inner product as a mapping 
TxT^C inherited from the inner product in H, as follows: 

oo oo N 

ffl*):=^f|$„ = £^C iAai V|*),|$)€T. (2.3) 

n— n— a — 1 

Here, ^ is the row-vector obtained by transposition and conjugation of the column-vector 
^> n , and ^l^n denotes the inner product of two elements of C N . By simplycity, we use 
the same notation (•]•) for the inner product in H and in T. In terms of this inner product, 
one can define a distance between |\ff) and |$) as usual: = - $|\& - < E>)] 1 / 2 , 

which converts T to a metric space. As this metric generates a (strong) topology in T, all the 
topological notions can be used for studying it, and it can be proved to be a complete metric 
space. From this point of view the right hand side in ( I2.2i is a Fourier series convergent in 
the strong topology in T. The fact that T is not a Hilbert space causes no difficulty from the 
physical point of view, since the underlying space H is a Hilbert space, and we introduce T 
only as a convenient tool. 

In the sequel we shall impose the following form for the action of the operator 7Yo 
introduced in equation ( 12. 1 1 on the basis \n, a): 

N 

H \n,a) = ^[D n>a!0 \n-l,f3)+D n+hajl3 \n+l,f3)+Q nta! i3\n,/3)U2A^ 
0=1 

This type of action corresponds to a tridiagonal matrix whose non-zero entries are A^- 
dimensional square matrices. The coefficients D n ^ a ^ and Q n ,a,/3 must satisfy D n>at p = 
D* nt f3 t u,Qn, a, p = Qn.p.a to ensure tne Hermitian character of Ho- 
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If for a fixed value of n we consider the N x N matrices with entries D n „ p and Qn, a, p 
as a matrix representation of the self-adjoint operators D n — 2?t and Q n = QT, whose action 
on is 

N N 

P»e a '■= Dn, a, &0 , Qn^a ■— Qn., a, &0 , (2.5) 

0=1 0=1 

then, according to ( 12. 4t . the action of the Hamiltonian 7i on the basis vector |n, a) takes the 
form of a three-term relation with operator- valued coefficients 

Ho\n, a) = T> n \n — 1, a) + V n+ i\n + l,a) + Q n \n, a) , T> \n, a) = 0, Vrt, Va. 

We shall allow the operator 7Yo to be unbounded in H, and to be defined on an appropriate 
dense set Bci, which is supposed to be invariant under the action of Hq. Then, V | S D 
the following Fourier series converges in H 

oo N 

H \^) = ^ X! l^n-l^n + Ipn+haDn+l + i>n,aQn]\n, a) . (2.6) 
n— a— 1 

With the help of this expression and < I2.2I > we can define the action of 7Y in T to be: 

oo 

Wo|4'> = V («o*)n |«>, 

(2-7) 

(Wo*)™ := P„+i*„+i +P„*„-i + Q„*„. 

Note that the right hand side in the first equation is a convergent Fourier series in T since it 
just corresponds to the convergent Fourier series J2.6I > in EL Moreover, Ho, as an operator 
acting in T, has a well-defined domain of definition D C T consisting of all elements of D 
considered as elements of T. We note that since D is invariant under the action of Ho, then D 
is also invariant under the action of Ho acting on T. 

The eigenvectors of Ho are obtained from the finite-difference eigenvalue problem 

+ Z>»*n-1 + Qn*„ = (2.8) 

with operator- valued coefficients (cf. (SJ). We do not indicate explicitly the dependence of 
\& n on £ in order to avoid a cumbersome notation. 

Let us suppose now that some solutions of the eigenvalue problem for Ho are known and 
we want to get a solution of the eigenvalue problem for another Hamiltonian Hi- According 
to the usual strategy of intertwining operators 1 9 1, one can try to solve this problem by finding 
an intertwiner (or transformation operator) C defined by 

CHo = HiC (2.9) 

where we consider all operators acting on the same domain D C T. Once C is found, the 
eigenkets |\&) of Hi are obtained by applying C to the eigenkets |\&) of Ho- I*) = ^\^)- 
Here, both |\&) and |^) correspond to the same eigenvalue £. 

We consider now a particular ansatz for C, by assuming that it acts on the basis element 
| n, a) of the space H as follows: 

N 

C\n,a) = ^[A ntai0 \n-l,[3) + B nia> p\n,[3)] (2.10) 
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with Aq a a = 0. We easily find that the action of C on any \^>) from U C H is 

oo 

£|*) = |*> = E* 



n=0 



n ,( 7 ; 

(2.11) 



= (£*)„ := A l+ i$ n+ i + 

Here A n and £>„ are operators acting on C N , which in the basis {e a } are represented by 
matrices with entries A n> Qi p and B n ^ Qi p. For completeness, we remark that Aq = 0. 

It should be mentioned that the second relations of ( 12. 71 and ( 12.1 1> define the action of 
Hq and L on the space of vector-sequences {^n}, which we shall denote by T£ 2 , since it 
is just a counterpart of the Hilbert space £ 2 . To be precise, we should remark that these are 
not really the same operators as the ones defined above, since they act on different spaces. 
Nevertheless, we will use the same notations for them and for the inner product in T£ 2 , 
believing that this will not cause any trouble to the reader. The inner product in T£ 2 is defined 
as usual by ( 12. 3t , where now := {^ n } , := {®n} are elements of T£ 2 . Once we have 
an inner product in T£ 2 , we are able to determine the adjoint of an operator C, defined as 
usually by the relation = $ e T£ 2 . We shall consider & as acting on 

finite elements which form a dense set in T£ 2 . Using equations ( I2.11> . we evaluate 
and find the action of £t 

(£t*)„ = ^t* n _i + fit* B (2.12) 
where A\ and B 7 \ are the operators adjoint to A n and B n , respectively. 

3. Discrete-matrix intertwining technique 

3.1. Intertwining technique in the space T£ 2 



The second relations of \2.1\ and ( 12. 1H define the action of TLa and C on any finite element 
{^ n } of T£ 2 , meaning that the sequence {^ n } has a finite length. For the product CHq one 
gets 

C(H *) n = A n+ i{V n+2 ^ n+2 + Q„+i*„+i +P„+i*„) 

+ B n {V n+l * n+1 + Q„V n + P„*„-i) . (3.1) 
A similar expression is obtained for HiC. Then the intertwining relation ( 12. 9> implies 

A n T> n+ i = VnAn+i (3.2) 

B n V n = V n B n ^ (3.3) 

An+l Qn+1 + B n V n+1 = V n+1 B n+1 + Q n An+l (3.4) 
A n+1 V n+1 + B n Q n = V n A n + QnB n . (3.5) 

When solving these equations (the unknowns being the operators X>„, Q n ,A n ,B n , and the 
data T> n , Q n ) we should take into account that Aq = T> o = 0. By assuming the operators 
A n , B n to be invertible, we easily find from (13.31 

V n = BMBn^y 1 (3.6) 

which being placed to d3.2t gives us 

A n V n+l {A n+l )- 1 = BnVniBn-x)- 1 . (3.7) 
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Similarly, from J3.5I > we have 

Q„ = An+iVn+^Bn)- 1 + B n Q„(S„) _1 - V n A n {B n )- x (3.8) 
which together with J3.4I > yields 

A n +\Q n +i + B n T> n+ i = T> n+ iB n +i + A n +iV n+ i(B n ) 1 A n +i (3.9) 

+ B n Q n {B n )- 1 A n+1 - V n A n (B n )- 1 A n+1 . 

To solve this equation we proceed to eliminate the variable B n by introducing a new auxiliary 
variable (actually, operator) a n defined by 

B n = An +1 a n . (3.10) 

In equation ( 13.91 we replace T> n by its expression ( 13. 61 , multiply the result from the left by 
(An+i) and use d3.7i and (13. lOi to rearrange the first term on the right hand side to get a 
nonlinear finite-difference equation where the only unknown is a n : 

Qn+i— r D n+2 (7 n+ i-'D n+1 a~ 1 — cr n (Q n ~'D n+ ia n —'D n a^ 1 : 1 )a~ 1 .(3.11) 

This equation having operator-valued coefficients may be interpreted as a matrix finite- 
difference counterpart of the derivative of a Riccati equation. It can be linearized and 
"integrated" by introducing the new variable U n defined by 

(Tn = -Un+lU- 1 . (3.12) 

Using it in ( 13. lit and simplifying, we get the difference equation G n +\ = G n , with 

G n = U' 1 [Qn + V n+1 U n+x U- 1 + V n U n -iU- x \U n , (3.13) 

whose solution is G n = A. The matrix A plays the role of an integration constant. This leads 
us to an equation for U n 

V n+l U n+1 + V n Un-l + QnU n = U n A (3. 14) 

which is identical to the eigenvalue problem ( 12. 8t . except for the fact that the vector-valued 
variable is replaced by the matrix- valued variable U n and the scalar parameter £ is 
replaced by the matrix A. We shall now show how solutions of (13. 14i may be obtained from 
solutions of d2.8l >. Let us choose A to be diagonal: A = diag(Ai, . . . Ajv). In fact, this is not a 
restriction, since this matrix is Hermitian (this follows from the Hermicity of the matrices T> 
and Q) and therefore it can be diagonalized by an appropriate unitary transformation. Denote 
the columns of the matrix U n by U\ tn , ■ ■ ■ , Um.u- If these vectors are solutions to the initial 
eigenvalue problem 

T> n+1 U i>n+1 + £>„Ui, n _i + Q n U hn = X l U hn , i = 1, . . . ,N (3.15) 

then the matrix U n is a solution to equation J3. 14b . Once the matrices U n are given, one can 
evaluate B n using (13. 121 and (13. 101 . provided the matrices A n are known. 

Now, we shall obtain an expression for A n . In order to do that, we have to solve equation 
i3.7\ , where B n are given in ( 13.1 Oi l, that is 

A n V n j r \A^ Srl = An+iRnA" 1 . (3.16) 

Note that the quantities 

R n = o n V n o~\, n = 1, 2, . . . (3.17) 

are known, because they depend only on the initial data T> n and the matrices a n (or U n ), 
which are given. Using very simple algebra, one obtains from ( I3.7I > 

V n+ iR n = (A n 1 A n +\R n ) 2 , R n V n+ i = {A n+1 A n T> n+ {f' . (3.18) 



Intertwining technique for a system of difference Schrddinger equations 



7 



It is a simple exercise to see that these equations are compatible, which means that equation 
( 13.161 1 has a nontrivial solution. Hence, using for instance the first one, we get 

An+i = A^Vn+M^R- 1 . (3.19) 

From this recursion relation we find all the matrices A n , n = 2, 3, . . . , in terms of Ai'. 



An — Al 



(V^) 1 / 2 ^ 1 } \{V s R 2 ) 1 ' 2 R^ 1 \---\(V n R n ^ 1 ) 1 ' 2 R- l _ 1 



Once A n and B n are found, from ( !3.21 -( l3~3T l we have two equivalent expressions for T> n : 

V n = An+lOnVnOn^A- 1 = A^n+l^l • (3.20) 

To find Q n we use ( 13.51 . which yields 

Qn = A n+1 (V n+1 + a n Q n ~ VnVnPn-l^-Klx ■ (3.21) 

It is necessary to mention that in finding the matrices V n and Q n we never required 
them to be self-adjoint. Therefore, our solution can provide us with both self-adjoint and non- 
self-adjoint matrices. The conditions T> n — T>\ and Q n = Q\ are really restrictions on the 
transformation function U n . 

Note that the quantities {U n } uniquely define both the transformation operator C and the 
new "potentials" T>„ and Q n . Therefore, we shall call U = {U n } the transformation function. 

3.2. A particular case 

In this Section we shall consider a particular nontrivial case corresponding to N = 2, where 
general formulas established in the previous Section take the simplest form. Moreover, 
we shall show that for this case one can establish one-to-one-correspondence between the 
spaces of solutions of the initial and transformed equations and investigate the factorization 
properties of transformation operators. 

Let us put N = 2 and choose the matrices T> n and Q n to be multiples of the 2x2 identity 
matrix X: T>„ = d n T, Q n = q n %- Let Uj n be real solutions to the scalar difference equation 

dn+lUj^L+l + dn"j,n-l + <7nWj,„ = XjUj >n , j = 1,2. (3.22) 

Then, the matrix-function 

U n =( Ul ' n ~ U2 ' n ) (3.23) 

is a solution to the initial eigenvalue problem ( I3.14t . with A = diag(Ai, A2). We will show 
that it is suitable to be used as the transformation function in our algorithm, and that it will 
lead us to a nontrivial exactly solvable matrix-difference eigenvalue problem. After finding 
the inverse function 

U-^A-'f U2 - n U2 - n ) . (3.24) 

where A„ = 2 ui tTl U2, n > we get from (13 . 1 21 

,-1 _ a -1 ' LOn 7r " 



a n = -U n+1 U-' = -A- (3.25) 

\ "n w n J 

with 

W„ = Ui tn+ iU2,n + Ul,nU2,n+l and 7T„ = Mi,„+iM 2 ,n - "l,nU2,n+l -(3.26) 

Two remarkable properties of the matrices (I3.25> should be mentioned: 
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(i) From equations J3.25l > and ( 13.261 1. and from the real character of Uj it is obvious that 
all the matrices a n are Hermitian: crj = a„; 

(ii) They form a commutative group with respect to matrix multiplication (this assertion can 
be trivially proved using ( I3.25» . 



The last fact simplifies considerably our formulas. For instance, all the matrices in the 
expression for A„ commute, and we easily find 

where A% has been chosen appropriately. The expression for B n follows from ( 13.101 

B n = -V 1 n / * 1 (U n+1 U- 1 ) 1 / 2 . 
Using J3.20i and I3.2H we deduce formulas for T> n and Q n : 

v n = (v n v n+1 Un-x u- 1 u n+1 u- 1 ) 1 ' 2 

Qn = Qn — Dn+l Un ^n+1 + "^n Mn-1 U n 1 . 

For our particular choice (T> n = d n I, Q n — q n Z), we get explicit expressions for these 



(3.27) 
(3.28) 

(3.29) 



matrices in terms of the known solutions of Q.221 . Indeed, using ( I3.23l l-l l3.24l i. and choosing 



the positive square root in the first of equations ( I3.29> we get T>„ = 



a+ 



where 



a± 



V d n d, 



■71+ 1 



!Ul, n +lU ln -l 



± 



I M2,7i+lW2,7i-l 



Similarly, from the second of equations J3.29I one finds Q„ = Q„ 



7 '-"71 

b+ = — 



U-l, 77,-1 i ^2,7t-l 

«l,n 



U2-. 



l n+l 



± 



6_ 

u 2 , 



Ml,7l+1 ^2,71+1 



(3.30) 



where 



(3.31) 



We notice that for U2 : n — ui^ n the matrices T> n and Q n become diagonal, with the nonzero 
elements coinciding with the previously obtained for the scalar case |5 1. 

From J3.30i and 13.3 II . it is clear that the transformed quantities are Hermitian provided 
a± and b± are real. On the other hand, it is easy to show that, under some reasonable 
assumptions on q n and d n , the operators A n and B n are anti-Hermitian. Indeed, from J3. 141 
we have 

V n+1 U n+x U- x + VMn-xU- 1 - U n KU~ l - Q n . (3.32) 

By choosing the matrix A appropriately, and provided the sequence {q n } is bounded from 
below, the matrices in both sides of d3 .321 may be negative definite for all values of n to get 
the matrix T> n+ ia n + T) n a n -i positive definite. This is achieved by taking the operators T> n , 
as well as a n , to be positive definite Vn, and the anti-Hermiticity of A n and B n follows from 
J3~27b and dl28t . 

Let us now calculate the superposition OH as an operator acting on Ti 2 . Using 
equations ( 12.1 1> . (12. 121 . J3.27I . J3.28I . and the anti-Hermiticity of A n and B n we find 

(£ f £*)„ = (V n+1 ^ n+1 + P„*„_i + Q n ^ n ) 

- (V n+1 U n+1 + V n Un-l + QnU n )U~ l ^ n . (3.33) 

The first term on the right hand side is simply the function (7^0^)71 (see equation ( 12. In the 
parentheses of the second term we see the action of the same operator on the transformation 
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function U n . Hence, using A3. 14i we transform the right hand side of the previous equation 
to the form = (Hq — UnMA^ 1 )^ n . For the particular case of the transformation 

function given in ( I3.23> we easily obtain 

U n MA-' = l( X : + ^ ^-^)=:A (3.34) 

and therefore 



2 \ Ai — A2 Ai + A2 



tf£ = H -A. (3.35) 

We observe here a difference between our Darboux technique and supersymmetric quantum 
mechanics usually based on the factorization C = Ho — A. The factorization constant (in 
our case A) does not coincide with an eigenvalue of the transformation function. We believe 
that this is a manifestation of the fact that supersymmetric quantum mechanics based only on 
the factorization of the Hamiltonian is a particular case of the Darboux transformation method, 
which in general allows to factorize other symmetry operators. A similar factorization takes 
place for the inverse superposition CC > = Hi — A which may be established by applying the 
operator C to formula d3.35i . 

Our results indicate the existence of a one-to-one correspondence between the solution 
spaces of the initial and transformed equations. To develop this, we find the second solution 
tl n °f d3.14> with a given value of A, T> n -\-iU n +i + V n U n -i + Q n U n — U n A. By eliminating 
Q n from this equation and using the adjoint form of A3. 14l >. we get U^T> n+ iU n+1 — 
U 7 ] +1 V n+1 U n + U^T> n U n -i — U^_ x T> n U n = U^U n A — AU^U n . Because of our particular 
choice of the transformation function, the right hand side of the last equation vanishes and the 
resulting difference equation can be "integrated" to yield 

ul_{D n U n - U^V n U n -i = U^V n+1 U n +i - U r l +1 V n+1 U n = W . 

From here we see that the matrix Wq plays the role of the Wronskian for our discrete 
eigenvalue problem, and we find the recursion relation forU n , which being iterated gives 

n 

U n = {Ul)- 1 U,t Uo+Y, V k' H)~ X ^ H-xT 1 Wo ■ (3.36) 
fe=i 

Now, we can act with C on this function in order to get a solution of the transformed equation 
with eigenvalue A. After some algebraic manipulation we arrive at 

S n EE (C U) n = V~l{ 2 (U n U^) 1 ' 2 (3.37) 

where the constant matrix Wo has been chosen to be the identity. We note also that the 
operator-valued function S n belongs to the kernel of the operator £>, i.e. C^S n = 0. This 
can be established easily by using formulas (12. 12i . J3.27i and J3.28i . Once S n is found, one 
can get the second solution S n of the transformed equation by using (13.361 with appropriate 
changes. The solution U n also corresponds to the matrix eigenvalue A. 

The operator C ' satisfies the intertwining relation adjoint to ( 12. 9> . This means that it 
is a transformation operator from solutions of the transformed equation to solutions of the 
initial one. Then, if E ^ Ai 2> the function = (C 1 ^) n is a nontrivial solution of the 
transformed equation, whereas $„ = (£^)„ is also a nontrivial solution, but to the initial 
equation. Solutions of the transformed equation with eigenvalues Ai,2 can be found using 
J3.36i — (I3.37i . Hence, we have constructed a one-to-one correspondence between the solution 
spaces of the initial and transformed eigenvalue problems. 
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4. Application 

We shall use the above developed discrete-matrix Darboux transformation to generate new 
exactly solvable two-channel potentials. Indeed, we will show that even starting with the 
simplest case of the two-channel uncoupled free particle Hamiltonian, we can obtain a family 
of nontrivial coupled exactly solvable potentials. In contradistinction to all known exactly 
solvable potentials, they are represented by infinite block-tridiagonal matrices. The eigenvalue 
problem for such matrices is not solvable by usual methods. 

Consider the free particle Hamiltonian ho = p 2 . Since the momentum operator p x is 
expressible in terms of the harmonic oscillator ladder operators a + = id/dx + ix/2 and 
a = id/dx — ix/2, ho = (a + a + ) 2 /4 is a quadratic form in a and a + . Therefore, the action 
of ho on the oscillator basis |n), which in coordinate representation is 

1> n (x) = (x\n) = (~t) n (nl2 n V2^)-^ 2 e- x2 / 4 H n (x/V2) (4.1) 
where H n (z) are Hermite polynomials, takes the form of a three term relation 

ho\n) = d n \n-2)+d n+2 \n+2)+q n \n), d n = V " ( 4 "~ 1) ; q n = f +£.(4.2) 

To derive J4.2i we have used a\n) = y/n\n — 1) and a + \n) = y/n + l\n + 1). Let \i/je) be 
a continuous spectrum eigenket of ho'. ho\ipE) = E\i/>e)- Then, using the Hermiticity of ho 
we get the scalar discrete eigenvalue problem 

+ d n+2 1p n +2 + Qni'n = Elj) n . (4.3) 

A "physical" solution ip n — tp n (E) — (iPe\ti),E > 0, to this equation can be easily obtained 
since it coincides with the Fourier image of the function (14. 1 1 : 

4> n {E) = 2(n\2 n V2^y 1/2 e- E H n {V2E) . (4.4) 

Consider now an uncoupled two-channel problem with the diagonal matrix Hamiltonian 
Ho = hoi- This is only the kinetic energy operator which acts in the Hilbert space 
H = H ®C 2 , with basis \n,a) = e a \n), e x = (1,0)*, e 2 = (0,1)*. The action of Ho 
on a two-component vector |W) = J2^Lo ^n\n), is given by (12. 7L where 

(Wo*)* - 2\+2*n+2 + ^„*„-2 + Qn^n £ C 2 (4.5) 

T^n = d n X and Q n = q n I. Since the solutions of the one-channel problem for ho ( 14. 3> are 
known, we know solutions of the two-channel problem for 7io- In particular, two-column 
vectors {^+{E)} = {(V>„(£), 0)*} and {^-{E)} = {(0, Vn(^))*} with V„(£) given in 
J4.4i are solutions to equation 

%+2*n+2 + ©„*n-2 + Qn*n = (4.6) 

with £ — E. The matrix U n given in ( I3.23> . with matrix elements 

u ln = (n\ 2 n )- 1 l 2 H n {^2\ 1 ) , u 2n = (n! 2 n )- 1 / 2 H n {^2\ 2 ) (4.7) 
A i. 2 < 0, is a solution to the matrix eigenvalue problem 

V n+2 U ri+2 + V n Un-2 + QnUn = U n K (4.8) 

with A = diag(Ai,A2). Remark that this two-channel problem is not exactly of the type 
considered in Section 2: here the eigenvalue problem involves shifting by two units the 
discrete variable n. Nevertheless, it is not difficult to see that we can consider here two 
independent problems, one for the even and the other for the odd values values of n. For the 
sake of simplicity we omit these details and give below only the final results. 
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(4.9) 



V\n, a) = Q n \n - 2, a) + G n +2\n + 2, a) + TZ n \n, a) , a = 1,2 

with coefficients (/„ and 7?.,,, that are to be determined. Let \^>) be an eigenvector of Hi, 
Hi\ty) — and { } be its Fourier coefficients over the basis \n). Let these quantities 
satisfy the similar discrete matrix eigenvalue problem 

Vj> n -2 + V n+2 V n+2 + Qn^n = . (4.10) 

Suppose also that the unspecified quantities Q n and lZ n are related to T> n and Q n by 

V n = g n + \jn(n- 1)1, Q n = H n + (q + \)T. (4.11) 

Consider the subclass of potentials d4.9t for which equation 14. 101 coincides with the Darboux 
transform of the initial eigenvalue problem J4.6i and express the transformation function U n 
in the form (13.231 with the entries given in 14.71 . All these assumptions are satisfied if the 
functions T> n and Q„ have the form 

(4.12) 

Qn = Qn - V n+2 U n Unl 2 + 2^4-2^ • 

The matrices (?„ and 7?. rl represent in this case the "potential differences" Q n = T> n — T> n , 
Tin = Qn — Qn which are expressed in terms of known quantities only. Solutions of 14. lOi 
are found with the aid of ( 12.1 H . A n and B n given in 13 .271 — (13 ,28> : 

*n = K+2^nU-l 2 ) 1/2 (^n+2 ~ Un+JG^n) • (4-13) 

Thus, we have solved the eigenvalue problem for a Hamiltonian Hi = Hq + V having in the 
basis \n, a) an interaction in the form of the block-tridiagonal matrix 

(k/3\V\na) = G fj , oi,n$n — 2,k 

(4.14) 

Gp, a ,n and Rp_ a>n are entries for Q n and lZ n respectively, and are found from (14.1 It — (H-.12I . 
Hence, using J3.30i and the notation H n (\/2Xi) = H n ,i and H n (y/2X~2) = i? n ,2, we find 



Qn = —\y/n(n -1)1 



0- 



a+ 



and lZ n = 



6_ 



with 



a± 



b± 



+ 2,1-^ti-2,1 I 



H n. + 2.2-H Vt-2,2 



(4.15) 



nin— 1) 
4 



H„- 



± 



(n+l)(rt+2) 
4 



± 



H„ 



H ns l Hn,2 

Note that the quantities H n i and -ff ri! 2 are real for n even and purely imaginary for n odd and 
therefore Q n and lZ n are real in either case. 

Let us find now the asymptotic behaviour of these quantities with respect to the discrete 
variable n. For this purpose we need an appropriate asymptotic expression for the Hermite 
polynomials, which can be found with the help of the relationship between Hermite and 
Laguerre polynomials, H 2n (z) = (-l)"2 2n n!L I T 1/2 (z 2 ), and using the series expansion for 
the Laguerre polynomials (see [10|,Eqs. 10.15(4), 10.15(5)) in which we keep only the two 
first terms and consider z 2 < 0, Im(z) > 0. We then have 



2"- 



-Inl 



1 



16n 



-3/2 >| 



(4.16) 
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and 

, — (-1\ n (In - TW! . „. . „,„i 

(4.17) 



H 2n (V2X) = ( 2) " ( ^ e*" 2 ^ [l - A + 0(n- 3 / 2 : 



From ( RTBl we obtain a+ = n/4 + 1/8 - l/(32n) + 0(l/n 2 ),6+ = 1/2 + 0(l/n 2 ), 
a_ = 0(l/n 2 ), and &_ = 0(l/n 2 ). The matrices Q n and 7?.„ become diagonal up to 
terms of order 1/n 2 . Moreover, up to quadratic terms in 1/n, a + coincides with d n+ \ and 
(<?n + &+) with <7n+i- This means that for n — > oo the Darboux transformed potentials X>„ 
and Q n coincide with the initial ones, up to the shifting n — > n + 1. This is the discrete 
multichannel analogue of a property of the usual Darboux transformation saying that usually 
the potential difference vanishes asymptotically. 

Now we will estimate the asymptotics of solutions of the transformed discrete eigenvalue 
problem. To do it we will assemble the column-vectors ^>^{E\) and ^>~(E2) in the matrix 
5 n , which is a matrix solution of the initial equation corresponding to the diagonal matrix 
eigenvalue with entries E\ and E 2 . The matrix 5„ := (^f^(Ei), ^~(E 2 )) (the columns 
(-^1,2) are given by J4. 1 3I > with the replacement of ty n by $^(£^2) respectively) is the 
corresponding matrix solution to the transformed equation 

H„ = VS n , V= V)H 2 {U n+2 U- 1 ) ^ 2 (W„W-^ 2 H„ +2 H- 1 - 1) . (4.18) 
Then, using the same asymptotics for the Hemite polynomials J4.17> . we obtain 



V P2 J q= |(V^- VA7) + o(n- 1 / 2 ). 

We see that for ri — > cxd the matrix "P does not depend on the variable n, a result which agrees 
with the fact that the potential difference vanishes asymptotically. This means, in particular, 
that if we choose the initial waves propagating only in one direction in both channels, the 
same asymptotic behavior will be exhibited by the solutions of the transformed equation. 
This property is precisely a two-channel analogue of the transparency for the transformed 
potential we have obtained. 



5. Conclusion 



In this paper we have introduced a difference intertwining operator for discrete Schrodinger 
equations with operator-valued coefficients and then we studied some of its basic properties, 
such as the factorization and the possibility to establish a one-to-one correspondence between 
the spaces of solutions of the initial and transformed equations. This approach allowed us to 
get new exactly solvable eigenvalue problems for block-tridiagonal matrices. By working out 
the particular example of a two channel Schrodinger equation, we have found a wide class of 
exactly solvable two-channel potentials represented by block-tridiagonal matrices. 

One of the possible applications of the method we have developed may be in describing 
the scattering of composite particles, such as nucleons or atoms, in the frame of the 
multichannel scattering theory |6|. This possibility is supported by the known applications 
of the J-matrix method for describing the one-channel scattering |3j|4). Precisely in this 
line, the next step of this work will be to find the transformation of the scattering data for the 
spectral problem on a semi-axis. 
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